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Abstract 


A  numerical  study  of  the  flow  of  liquid  propellant  (LP)  in  a  regenerative  LP  gun  (RLPG)  is 
described.  The  model  simulates  the  flow  of  LP  from  the  liquid  reservoir,  through  the  orifice,  and 
into  the  combustion  chamber.  The  model  is  based  on  a  space-time  finite  element  method  and 
can  automatically  handle  the  deformation  of  the  computational  domain.  A  mesh-moving  scheme 
is  used  to  update  the  mesh  at  every  time  step. 

Two  shots  of  the  30-mm  RLPG  are  studied.  One  is  a  small-charge  shot,  and  the  other  is  a 
medium-charge  shot.  The  results  from  the  simulations  compare  well  with  data  from  another 
numerical  model  of  the  RLPG  and  with  experimental  data.  The  simulations  provide  qualitative 
details  of  the  transient  phenomena  that  occur  in  the  orifice  during  the  firing  cycle.  The 
agreement  of  the  results  from  the  model  with  experimental  data  provide  confidence  in  the 
accuracy  of  the  model.  The  numerical  model  will  be  used  to  study  several  RLPG  shots  under 
a  variety  of  conditions. 
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1.  Introduction 


The  regenerative  liquid  propellant  gun  (RLPG)  is  being  studied  by  the  U.S.  Army  as  an 
advanced  technology  for  future  gun  systems.  There  are  several  technical  issues  that  must 
be  resolved  before  RLPG  technology  is  sufficiently  mature  for  a  fielded  weapon  system.  In 
order  to  help  researchers  understand  the  gun’s  firing  cycle  sufficiently  to  solve  the  technical 
challenges  presented  by  the  RLPG,  30-mm,  105-mm  and  155-mm  RLPGs  have  been  built 
and  tested.  In  addition,  a  great  deal  of  numerical  modeling  of  the  firing  cycle  has  been 
performed  in  support  of  this  goal. 

The  firing  cycle  of  the  concept  VIC  (6-C)  RLPG  begins  with  the  liquid  propellant  (LP) 
stored  in  a  liquid  reservoir  (see  Figure  1).  There  are  two  pistons:  (1)  the  inner,  or  control, 
piston,  which  moves  along  the  centerline  of  the  gun;  and  (2)  the  outer,  or  injection,  piston, 
which  surrounds  the  inner  piston.  An  external  ignitor  injects  hot  gas  into  the  combustion 
chamber  to  start  the  firing  cycle.  As  the  pressure  in  the  chamber  increases,  the  pistons  are 
driven  toward  the  breech.  There  is  a  damper  that  sits  behind  the  control  piston  and  slows 
the  piston  near  the  end  of  stroke.  Early  in  the  firing  cycle,  the  damper  applies  relatively 
little  force  to  the  control  piston,  while,  at  the  same  time,  the  motion  of  the  injection  piston 
is  slowed  by  the  presence  of  LP  in  the  reservoir.  Therefore,  the  control  piston  moves  a  little 
faster,  early  in  the  firing  cycle,  than  the  injection  piston,  and  the  difference  in  the  motion  of 
the  pistons  opens  an  annular  orifice  between  the  pistons. 

Throughout  most  of  the  firing  cycle,  the  ratio  between  the  pressure  in  the  reservoir 
and  the  pressure  in  the  chamber  is  equal  to  the  inverse  of  the  ratio  between  the  effective 
areas  over  which  the  pressures  act  on  the  injection  piston.  In  other  words,  the  force  that 
the  chamber  pressure  exerts  on  the  injection  piston  is  roughly  equal  to  the  force  that  the 
reservoir  pressure  exerts  on  the  piston.  Since  the  ratio  of  the  area  on  the  chamber  side  to 
the  area  on  the  reservoir  side  is  1.47  for  the  30-mm  RLPG,  the  pressure  in  the  reservoir  is 
roughly  1.47  times  the  pressure  in  the  chamber  during  most  of  the  firing  cycle.  The  higher 
pressure  in  the  reservoir  drives  the  LP  from  the  reservoir  through  the  orifice  and  into  the 
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Figure  1.  Sketch  of  the  Concept  VIC  RLPG. 


combustion  chamber.  The  LP  combusts  in  the  chamber,  accelerating  the  projectile  down 
the  gun  tube. 

Occasionally,  the  firing  cycle  does  not  proceed  properly.  The  injection  piston  will  slow 
or  reverse  its  direction  of  travel  before  the  end  of  stroke.  This  is  called  a  piston  hesitation 
or  reversal,  or  simply  a  reversal.  Not  all  causes  of  reversals  have  been  identified.  One  cause 
may  be  combustion  in  the  orifice  or  reservoir,  and  reversals  of  this  type  are  the  focus  of  this 
work.  The  source  of  the  combustion  in  the  reservoir  is  not  known,  although  one  possibility 
is  a  disruption  of  the  flow  in  the  orifice. 

In  order  to  test  this  hypothesis,  a  better  understanding  of  the  flow  of  LP  in  the  orifice 
is  needed.  For  instance,  it  is  not  yet  known  if  the  flow  separates  from  the  surface  of  the 
pistons,  and  what  conditions  cause  this  flow  separation.  As  a  part  of  the  study  of  the  RLPG 
firing  cycle,  a  30-mm  RLPG  was  built  and  studied  [1]  and  is  still  used  as  a  tool  in  the  study 
of  the  RLPG  firing  cycle.  While  the  30-mm  RLPG  is  a  very  valuable  experimental  fixture, 
it  is  prohibitively  difficult  to  measure  enough  experimental  data  in  the  orifice  to  describe  the 
flow  phenomena  that  occur  there. 

With  the  limitations  that  exist  on  the  experimental  study  of  the  gun,  it  is  necessary  that 
numerical  modeling  provide  at  least  a  qualitative  description  of  the  flow  phenomena  that 
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occur  during  the  firing  cycle.  Numerical  modeling  has  been  used  to  study  the  firing  cycle 
of  the  RLPG  for  a  number  of  years.  Work  by  Coffee  et  al.  [1-6]  has  been  focused  primarily 
in  two  areas.  The  first  is  the  lumped  parameter  modeling  of  the  RLPG  firing  cycle.  The 
second  involves  modeling  the  combustion  chamber. 

In  the  lumped  parameter  approach,  the  spatial  variations  of  the  pressure,  density,  and 
temperature  within  both  the  reservoir  and  the  chamber  are  neglected.  Also,  the  flow  through 
the  orifice  is  computed  using  an  assumption  of  Bernoulli  flow  [1].  While  the  data  from  the 
lumped  parameter  model  match  the  experimental  data  very  well  and  are  very  valuable  to 
researchers,  they  do  not  provide  many  more  details  of  the  flow  phenomena  in  the  orifice  than 
the  experimental  data  provide. 

The  combustion  chamber  model  [6]  simulates  the  breakup  of  the  jet  of  LP  entering  the 
chamber  and  the  combustion  of  the  resulting  droplets.  It  captures  the  pressure  waves  as 
they  move  around  the  chamber,  causing  pressure  oscillations  at  the  chamber  end  of  the 
orifice.  The  mass  flow  rate  of  the  LP  into  the  chamber,  however,  is  computed  by  the  lumped 
parameter  model,  and  the  shape  of  the  velocity  profile  across  the  orifice  is  approximated. 

The  current  work  involves  the  modeling  of  the  liquid  reservoir  and  the  orifice  in  the 
30-mm  RLPG.  The  goal  of  this  work  is  the  development  of  a  model  that  can  capture  the 
features  of  the  flow  in  the  orifice.  The  model  used  in  this  work  was  introduced  previously  [7]. 
Because  of  the  complicated  shape  of  the  RLPG  interior,  the  use  of  unstructured  meshes  and 
an  automatic  mesh  generator  makes  discretization  of  the  computational  domain  much  easier. 
The  finite  element  method  was  chosen  for  the  model  because  of  its  ability  to  automatically 
handle  unstructured  meshes. 

Also,  the  motion  of  the  pistons  changes  the  size  and  shape  of  the  domain;  therefore, 
the  Deformable-Spatial-Domain/Stabilized-Space-Time  (DSD/SST)  formulation  [8-10]  of 
the  governing  equations  is  used.  In  the  DSD/SST  formulation,  a  finite  element  discretization 
of  both  space  and  time  is  performed.  Because  of  this,  the  changes  in  the  size  and  shape  of 
the  computational  domain  can  be  handled  automatically,  even  if,  as  in  this  case,  the  changes 
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are  not  know  ahead  of  time  and  must  be  computed  as  part  of  the  overall  solution.  The 
DSD/SST  formulation  works  by  treating  time  as  another  spatial  dimension.  The  cost  is 
that  a  two-dimensional  spatial  problem  becomes  a  three-dimensional  space-time  problem. 
To  minimise  this  cost,  the  space-time  domain  is  subdivided  into  space-time  slabs,  analogous 
to  time  steps  in  a  semi-discrete  formulation.  The  DSD/SST  formulation  is  then  solved  one 
space-time  slab  at  a  time. 

Two  stabilization  methods  are  used  in  the  simulation  of  the  RLPG.  The  first  is  the 
streamline-upwind/Petrov-Galerkin  (SUPG)  method  [11],  which  controls  the  instability 
caused  by  advection-dominated  flows.  It  has  been  used  in  modeling  a  variety  of  compressible 
fluid  flow  problems  [12,13].  The  second  stabilization  method  is  a  discontinuity  capturing 
method  [14],  which  prevents  overshoot  and  undershoot  in  the  neighborhood  of  shocks  and 
sharp  gradients.  It  also  has  been  used  in  a  variety  of  flow  problems  [13, 15].  While  there  are 
no  shocks  in  the  interior  of  the  RLPG,  there  are  sharp  shear  layers  between  the  jet  of  LP 
entering  the  combustion  chamber  and  the  rest  of  the  LP  in  the  chamber. 

2.  Governing  Equations 

In  the  numerical  model  discussed  in  this  paper,  the  propellant  is  assumed  to  be  single 
phase  throughout  the  gun.  The  propellant,  which  in  the  RLPG  is  XM46,  must  be  assumed 
to  be  compressible,  because  of  the  high  pressure  in  the  gun.  It  is  often  convenient  to  use  the 
bulk  modulus  to  compute  the  compressibility  of  a  liquid.  The  bulk  modulus  is  defined  as 

(i) 

and  is  assumed  to  be  a  linear  function  of  the  pressure,  i.e., 


I<  =  Ki  +  Ktp. 


(2) 


In  these  equations,  K  is  the  bulk  modulus  and  p,  p ,  and  8  are  the  density,  pressure  and 
temperature  of  the  liquid.  Also,  K\  and  Ki  are  experimentally  determined  constants.  For 
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any  situation  in  which  the  pressure  of  the  liquid  is  within  an  order  of  magnitude  of  the  bulk 
modulus  (as  is  the  case  in  the  RLPG),  the  compressibility  of  the  liquid  cannot  be  neglected. 


Equations  (1)  and  (2)  can  be  combined  to  generate  the  equation  of  state  for  this 
compressible  liquid, 


(3) 


The  density  at  gage  pressure  of  zero  is  po-  The  pressure  of  the  liquid  is  a  function  of  only 
density;  therefore,  the  liquid  is  said  to  be  barotropic. 


The  dynamics  of  the  liquid  in  the  model  are  governed  by  the  compressible  Navier-Stokes 
equations.  Since  the  propellant  is  barotropic,  and  its  viscosity  is  assumed  to  be  constant, 
the  conservation  of  energy  equation  decouples  from  the  conservation  of  mass  and  momentum 
equations.  Only  the  latter  are  required  to  describe  the  flow  of  the  liquid.  The  former  can 
be  solved  using  data  from  the  other  equations,  if  desired.  The  conservation  of  mass  and 
momentum  equations  can  be  written  in  conservation  law  form  as  Vx  €  Qt,t  €  (0,  T) 


l  +  V-tpu)  =  0 


(4) 


d(pu) 

dt 


+  V  •  (upu)  +  Vp 


VT, 


(5) 


where  t  is  time,  (0,  T )  is  the  time  span  of  interest,  and  f lt  is  the  spatial  domain.  The  subscript 
“t”  implies  that  the  spatial  domain  may  be  changing  with  time.  Tt  is  the  boundary  of  Qt. 
Also,  u  is  the  fluid  velocity,  and  T  is  the  Newtonian  viscous  stress  tensor,  defined  as 

T  =  ,,  (Vu  +  (Vuf  -  j(V-u)l),  (6) 

where  p  is  the  viscosity  of  the  liquid,  and  I  is  the  identity  tensor. 
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The  geometry  of  the  RLPG  and  the  flow  features  of  interest  are  axisymmetric.  Assuming 
axisymmetry  of  the  solution,  the  governing  equations,  equations  (4)  and  (5),  are  written  in 
conservation  variables  form  as  Vx  €  €  (0,  T) 

The  radial  and  axial  coordinates  are  denoted  r  and  z.  U,  Fz,  and  Fr  are  the  vector  of 
conservation  variables  and  the  Euler  flux  vectors,  respectively.  They  are  defined  as 


U  =  puz  ,  Fz  =  uzpuz  +  p  ,  and  Fr  —  uTpuz 

,  pUr  /  l  UzpUr  /  l  UrpUr  +  p 


Sr  =  — U. 
r 


E  contains  the  dissipative  terms  and  is  defined  as 


Equation  (7)  can  be  written  in  quasi-linear  form  as:  Vx  G  f lt,i  €  (0,T) 


dU  A  dJJ  t  d\J  0 

~k~  +  — f-  Ar-^  K  Sr  —  V  •  E  —  0. 

at  oz  or 


Az  and  Ar  are  the  Jacobians  of  Fz  and  Fr  with  respect  to  U, 


A  _dFz  _8Fr 
z  dv'  T  dir 


Appropriate  boundary  and  initial  conditions  are  added  to  equation  (11)  to  complete  the 
problem. 
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3.  Numerical  Model 


The  initial  geometry  of  the  RLPG  interior  is  fairly  complex.  In  addition,  because  the 
pistons  move  a  large  amount  during  the  firing  cycle,  a  new  mesh  must  be  generated  from  time 
to  time  during  the  simulation  to  prevent  excessive  mesh  distortion.  For  these  reasons,  an 
automatic  mesh  generator,  which  uses  Delaunay- Voronoi  methods  [16],  is  used  to  discretize 
the  spatial  domain  using  an  unstructured  mesh,  which  consists  of  three-node  triangular 
elements.  These  triangular  elements  in  space  are  made  into  six-node  prism  elements  in  the 
space-time  domain  by  adding  the  temporal  dimension  to  each  element. 

The  space-time  domain  is  subdivided  into  space-time  slabs  by  partitioning  the  temporal 
dimension  (0,T)  into  subintervals  /„  =  (tn,tn+i),  where  0  =  t0  <  ti  <  •  •  •  <  ts  =  T.  The 
time  slab  corresponding  to  In  (denoted  Qn)  is  enclosed  by  Dn,  Dn+1 ,  and  Pn  (where  Dn  =  Dtn, 
and  Pn  is  the  surface  defined  by  T*  as  t  varies  from  tn  to  tn+ 1).  Also,  Pn  is  decomposed 
into  (P„)g  and  (P„)h,  the  parts  of  the  boundary  on  which  axe  imposed  Dirichlet-type  and 
Neuman-type  boundary  conditions,  respectively. 

The  function  spaces  for  the  trial  solutions  and  the  weighting  functions,  and  V£,  are 
defined  on  each  space-time  slab.  The  functions  are  piecewise  linear  and  continuous  in  space 
and  piecewise  linear  in  time,  but  discontinuous  across  the  interface  between  neighboring 
space-time  slabs.  The  discontinuity  of  the  shape  functions  in  time  means  that  the  mesh  used 
in  one  space-time  slab  need  not  be  used  in  the  next  slab  (i.e.,  remeshing  can  occur  across 
the  slab  interface). 
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The  DSD/SST  formulation  of  equation  (11)  is,  given  (Ufc)“,  find  UA  €  S*,  such  that 


VWA  € 


+  /  VWA  :  E hdQ  +  [  (W*)+  •  ((Ufc)+  -  (U*):) 

jQn 


dn 


dWh\  fdVh  .  AhdVh  ,  ihdUh 


+  Ai^  +  A*^-  +  sMdg 


(n«i)»  ,  (  f)Wh  ,  ^<9W 

+  S  LTh  -^  +  {Ai)  ^ 

(ry  f  gf  dW^  dQ  =  /  W'1  H'1  dP.  (13) 

Si  V  V  3*  dz  +  dr  dr  ;  v  V.)H 


In  this  formulation, 


(U‘)J  =  limU(in±£), 

[  ( ...)dQ  =  f  / 

JQn  Jin  J&n 

[ J...)JP  =  [  J  (...)dTdt. 

J  Pn  J  In  J^n 


(14) 

(15) 

(16) 


Equation  (13)  is  applied  sequentially  to  Qo,  Q i,  •  •  • ,  Qn- i,  with 


(uk)0-  =  us, 


(17) 


where  Uo  is  the  initial  value  of  U. 

The  time-discontinuous  Galerkin  formulation  of  equation  (11)  consists  of  the  first  three 
terms  on  the  left-hand  side  and  the  term  on  the  right  of  equation  (13).  The  third  term  on 
the  left  weakly  enforces  continuity  of  the  solution  across  slab  interfaces.  The  SUPG  term  is 
the  fourth  term  on  the  left,  consisting  of  a  series  of  element  level  integrals.  The  discontinuity 
capturing  term  is  the  fifth  term  on  the  left. 
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The  discontinuity  capturing  coefficient,  S ,  is  [17] 


1 

2 

(18) 


The  SUPG  coefficient,  r,  is  a  diagonal  matrix  defined  as  [10] 


6  = 


a^  +  a;^+s;  ii? 


'A'1 

Ao 


Ensd 

j= 1 


dU 

dxz  dxz 


2~  + 
A-1  ^ 


96  9U 
dXr  dxr 


2_ 

Ar1 


r  =  max[0,  r0  -  ts  -  Td ]. 


(19) 


The  contribution  from  the  advection  terms  is  ra;  Ts  and  Td  are  corrections  for  the 
discontinuity  capturing  term  and  the  physical  diffusion.  Without  the  correction  terms,  the 
solution  is  overly  diffuse.  The  terms  in  equation  (19)  are  defined  as 


_  /tmax(  1/3,|) 
2(c+  |u-/?|)  ’ 

=  _ £ _ i 

(C+|U-/?|)2  ’ 

(c+  |u-/3|)»  ' 


(20) 

(21) 

(22) 


The  element  length  and  the  speed  of  sound  in  the  fluid  are  denoted  h  and  c.  Also, 


and 


V II  U  \\l 

ii  v  ii  u  in. 


o 

o 

o 

o 

o 

o 

II 

o 

o 

and  Kr  = 

o 

o 

o 

o 

^  0  0  &  J 

(23) 


(24) 
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In  equation  (13), 


VW*  :  E* 


dwrduT 

2/,arar 


2  dwr 
3^ 


(V-u) 


+t* 


.  <9u;*  duz  2  . 

+2"^71j7-3'‘^7(v'u) 

+2/1——  -  |/i— (V  ■  u), 


r  r 


(25) 


where 


V  •  u  = 


dur  ur 

17+  s7  +  7’ 


(26) 


and  and  uv  denote  the  weighting  functions  corresponding  to  the  axial  and  radial 
momentum  balance  equations. 

During  the  simulation  of  the  firing  cycle,  the  mesh  moves,  based  on  the  motion  of  the 
pistons.  An  automatic  mesh-moving  scheme  [18]  treats  the  computational  domain  as  an 
elastic  solid.  The  motion  of  the  pistons  is  used  as  input  for  the  mesh-moving  scheme,  which 
then  computes  the  new  location  of  the  interior  nodes.  In  order  to  prevent  excessive  distortion 
of  the  mesh,  a  new  mesh  is  generated  periodically  using  the  automatic  mesh  generator. 

Equation  (13)  generates  a  system  of  coupled  nonlinear  equations  in  each  time  slab.  The 
equation  system  is  solved  using  the  Newton- Raphson  method.  The  Generalized  Minimal 
Residual  (GMRES)  algorithm  [19],  with  a  diagonal  preconditioner,  is  used  to  solve  the  linear 
equation  system  generated  by  the  Newton-Raphson  method. 

When  the  code  enters  a  space-time  slab,  the  flow  variables  in  the  new  space-time  slab 
are  set  equal  to  the  flow  variables  in  the  previous  slab.  The  motion  of  the  injection  piston  is 
estimated  using  the  fluid  data  from  the  previous  time  slab,  and  the  mesh  is  updated,  based 
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on  the  piston  motion.  The  fluid  simulation  is  updated  using  the  new  mesh  and  piston  motion. 
A  second  iteration  is  performed  between  the  piston  motion,  the  mesh,  and  the  flow  variables. 
In  general,  two  iterations  between  the  three  parts  of  the  code  within  each  space-time  slab 
provide  sufficient  convergence. 

4.  Model  of  LP  Flow  in  the  30- mm  RLPG 

The  DSD/SST  formulation  (equation  (13))  was  used  to  model  the  flow  of  the  LP  during 
test  firings  of  the  30-mm  RLPG.  At  the  beginning  of  the  firing  cycle,  the  pistons  are  in 
contact,  separating  the  reservoir  from  the  combustion  chamber.  The  ability  to  handle  the 
merging  of  two  computational  domains  into  one  as  a  part  of  the  simulation  has  not  been 
added  to  this  model.  Therefore,  the  pistons  are  slightly  separated  at  the  beginning  of  the 
simulation.  The  opening  is  so  small,  however,  that  only  one  element  spans  the  gap.  Since 
there  are  no-slip  boundary  conditions  on  the  piston  surfaces,  no  LP  can  flow  through  the 
orifice,  as  though  the  pistons  were  in  contact.  This  approximation  has  worked  quite  well. 

As  mentioned  previously,  the  LP  is  assumed  to  be  a  single  phase  throughout  the  LP 
interior;  no  combustion  is  included  in  the  current  model.  Since  the  combustion  chamber  in 
the  actual  gun  is  primarily  filled  with  gas  resulting  from  combustion,  the  current  model  does 
not  properly  simulate  the  flow  in  the  chamber.  Therefore,  while  the  computational  domain 
includes  the  entire  liquid  reservoir  and  orifice,  it  does  not  include  the  entire  combustion 
chamber  (see  Figure  2).  The  boundary  that  cuts  the  combustion  chamber  moves  with  the 
pistons  during  the  firing  cycle.  To  include  the  effect  of  combustion  in  the  simulation,  the 
pressure  in  the  combustion  chamber  is  taken  from  experimental  data  or  a  lumped  parameter 
simulation  of  the  shot  and  imposed  on  the  boundary  of  the  computational  domain  in 
combustion  chamber. 

The  motion  of  the  injection  (outer)  piston  is  computed,  based  on  the  pressure  forces 
from  the  LP  acting  on  it.  Viscous  and  friction  forces  on  this  piston  are  neglected.  The 
injection  piston  is  not  allowed  to  overtake  the  control  piston.  The  damper,  which  sits  behind 
the  control  piston,  is  not  modeled.  To  include  the  effect  of  the  damper,  the  motion  of  the 
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Boundary  of 

Computational 
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Figure  2.  The  Computational  Domain  Used  in  This  Study. 

control  piston  is  taken  from  experimental  data  or  a  lumped  parameter  simulations  of  the 
shot.  Both  of  the  pistons  act  as  rigid  bodies  and  only  move  axially. 

The  model  of  the  RLPG  assumes  axisymmetry.  A  no-slip  boundary  condition  is  imposed 
on  the  solid  surfaces,  and  a  no-penetration  boundary  condition  is  imposed  on  the  symmetry 
axis.  The  reservoir  is  initially  prepressurized,  while  the  chamber  is  initially  at  ambient 
pressure.  The  initial  condition  of  the  pressure  in  the  simulation  reflects  the  prepressurization. 
The  initial  pressure  varies  linearly  along  the  length  of  the  orifice  to  prevent  a  stability  problem 
caused  by  having  a  jump  in  the  initial  pressure  at  one  point. 

5.  30- mm  RLPG  Test  Firings 

The  30-mm  RLPG  was  fired  at  the  U.S.  Army  Research  Laboratory  at  Aberdeen 
Proving  Ground,  MD.  The  LP  used  was  XM46,  a  monopropellant  composed  of  hydroxyl 
ammonium  nitrate  (HAN),  triethanol  ammonium  nitrate  (TEAN)  and  water.  The  constants 
in  equation  (3)  for  XM46  are  K\  =  5350  MPa,  K2  =  9.11,  and  po  =  1430  Two 
shots  of  the  30-mm  RLPG  were  simulated:  shots  38  and  128. 

5.1  Shot  38.  Shot  38  was  fired  using  94.3  cm3  of  LP  and  an  initial  chamber  volume  of 
290  cm3  [20].  While  the  chamber  and  reservoir  pressure  data  from  this  shot  are  considered 
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to  be  accurate,  the  piston  travel  data  are  suspect.  Therefore,  the  data  used  as  input  for  this 
simulation  came  from  a  lumped  parameter  model  of  the  shot,  which  agreed  quite  well  with 
the  experimental  pressure.  The  mesh  of  the  initial  computational  domain  and  a  closeup  of 
the  mesh  in  the  orifice  are  shown  in  Figure  3.  The  initial  mesh  consists  of  3,029  spatial 
nodes  and  5,557  elements.  As  the  firing  cycle  proceeds  and  the  pistons  move,  a  new  mesh 
is  periodically  generated.  In  this  simulation,  a  new  mesh  was  generated  every  10  time  steps 
(the  time  step  size  was  0.0167  ms)  until  12.6  ms,  when  a  new  mesh  was  generated  every  time 
step  because  the  injection  piston  and  the  back  wall  of  the  reservoir  were  so  close  together. 


Figure  3.  Shot  38:  Mesh  at  Time  4.4  ms. 
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Time  t  =  0  was  set  to  be  the  beginning  of  the  ignition  sequence,  which  lasted  until  4.4  ms. 
The  simulation  started  once  the  ignition  sequence  ended,  and,  in  keeping  with  the  timing  of 
the  experiment,  the  initial  time  of  the  simulation  was  set  to  4.4  ms.  The  meshes  at  times 

6.4  ms,  8.4  ms,  10.4  ms,  and  12.4  ms  are  shown  in  Figures  4  and  5.  The  meshes  in  the  orifice 
are  shown  in  Figure  6.  The  meshes  used  during  the  simulation  of  shot  38  ranged  in  size 
from  5,600-14,500  elements  and  3,000-7,500  spatial  nodes.  The  meshes  shown  in  Figures  4 
and  5  have  5,176,  5,445,  5,575,  and  4,526  spatial  nodes  and  9,851,  10,389,  10,649,  and  8,591 
elements,  respectively. 

The  computational  domain  is  noticeably  larger  at  8.4  ms  than  at  10.4  ms,  as  can  be  seen 
in  Figures  4  and  5,  but  the  meshes  have  roughly  the  same  number  of  nodes  and  elements.  As 
can  be  seen  in  Figure  6,  the  orifice  at  both  times  is  nearly  the  same  size.  Since  the  majority 
of  the  nodes  and  elements  are  in  the  orifice,  there  axe  about  as  many  nodes  and  elements  at 

8.4  ms  as  at  10.4  ms. 

Control  piston  motion  and  chamber  pressure  data  from  a  lumped  parameter  simulation 
of  shot  38  were  used  as  input  for  the  DSD/SST  simulation.  In  the  test  firing,  the  reservoir 
pressure  was  measured  in  the  middle  of  the  back  wall  of  the  reservoir.  Therefore,  the  pressure 
predicted  by  the  model  for  the  middle  of  the  back  wall  of  the  reservoir  is  recorded  and 
compared  to  the  lumped  parameter  data,  which  agree  well  with  the  experimental  data.  The 
reservoir  pressure  from  experimental  data,  the  current  model,  and  the  lumped  parameter 
model,  and  the  chamber  pressure  from  the  experimental  data  and  the  lumped  parameter 
model  are  shown  in  Figure  7. 

The  agreement  between  the  results  of  the  DSD/SST  model  and  the  lumped  parameter 
model  is  good,  and  both  models  agree  quite  well  with  the  experiment.  The  oscillations  early 
in  the  firing  cycle  are  caused  by  the  piston  and  the  column  of  fluid  in  the  reservoir  acting  as 
a  mass  and  spring.  The  oscillations  in  the  data  from  the  DSD/SST  model  late  in  the  firing 
cycle  are  caused  by  remeshing.  This  is  corrected  in  the  later  simulations  by  using  a  smaller 
time  step.  A  time  step  size  of  0.01  ms  is  sufficient,  resulting  in  a  new  mesh  being  generated 
every  0.1  ms,  instead  of  every  0.167  ms. 
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Figure  4.  Shot  38:  Global  View  of  Mesh  at  Times  6.4  ms  and  8.4  ms. 

The  injection  piston  motion  from  the  DSD/SST  and  the  lumped  parameter  simulations  of 
shot  38  and  the  experimental  data  and  the  control  piston  motion  from  the  lumped  parameter 
simulation  and  the  experimental  data  are  shown  in  Figure  8.  The  predicted  piston  travels 
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Figure  5.  Shot  38:  Global  View  of  Mesh  at  Times  10.4  ms  and  12.4  ms. 


do  not  agree  well  with  the  experimental  data.  The  injection  piston  motion  data  from  the 
numerical  models  agree  reasonably  well,  although  some  differences  are  seen.  The  lumped 
parameter  model  assumes  that  flow  through  the  orifice  is  developed  immediately,  while  in  the 
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Figure  6.  Shot  38:  Mesh  in  the  Orifice  at  Times  6.4  ms,  8.4  ms,  10.4  ms,  and 
12.4  ms. 


DSD/SST  model,  the  flow  needs  time  to  develop.  This  causes  the  early  difference  in  injection 
piston  travel  seen  in  Figure  8;  in  the  DSD/SST  model,  the  injection  piston  is  slowed  while 
flow  m  the  orifice  develops.  Later  in  the  firing  cycle,  the  injection  piston  in  the  DSD /SST 
model  nearly  catches  up  with  the  injection  piston  in  the  lumped  parameter  simulation. 


17 


4  5  6  7  8  9  10  11  12  13  14 

Time  (ms) 

Figure  7.  Shot  38:  Pressure  at  the  Back  Wall  of  the  Reservoir: 

DSD/SST  Prediction  (Line),  Lumped  Parameter  Model  (Dash),  and 
Experimental  Data  (Long  Dash).  Pressure  in  the  Combustion 
Chamber:  Lumped  Parameter  Model  (Dot)  and  Experimental 
Data  (Dot-Dash). 

The  last  test  of  the  accuracy  of  the  DSD/SST  model  is  the  discharge  coefficient  of  the 
orifice,  defined  as 


CD 


(Pr  ~  !>7) 


(27) 


The  mass  flow  rate  of  the  LP  through  the  orifice  is  denoted  my;  pr,  pc,  and  pr  are  the 
pressure  in  the  reservoir  and  the  chamber  and  the  density  in  the  reservoir,  respectively. 
A  is  the  cross-sectional  area  of  the  orifice.  The  generally  accepted  value  for  the  discharge 
coefficient  is  0.95  [1].  The  discharge  coefficient  is  computed  by  the  DSD/SST  model  and  is 
shown  in  Figure  9.  During  the  first  2  ms  of  the  simulation,  the  flow  is  developing  and  the 
discharge  coefficient  rises  rapidly.  It  levels  off  between  0.85  and  0.9,  which  is  a  little  lower 
than  expected.  This  means  that  the  mass  flow  rate  of  LP  per  unit  area  through  the  orifice 
computed  by  the  model  is  likely  a  little  lower  than  what  actually  occured  in  shot  38. 
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Figure  8.  Shot  38:  Motion  of  the  Injection  Piston:  DSD/SST  Prediction  (Line), 
Lumped  Parameter  Model  (Dash),  and  Experiment  (Long  Dash). 
Motion  of  the  Control  Piston:  Lumped  Parameter  Model  (Dot)  and 
Experiment  (Dot-Dash). 
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Figure  9.  Shot  38:  Discharge  Coefficient  in  the  Orifice:  DSD/SST  Prediction. 


The  low  discharge  coefficient  means  that  the  model  is  predicting  more  kinetic  energy 
loss  in  the  orifice  than  occurs  in  the  physical  gun.  The  method  of  the  kinetic  energy  loss  is 
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probably  its  conversion  to  heat  by  viscous  effects.  In  the  model  used  in  this  work,  the  drivers 
of  the  problem  are  the  motion  of  the  control  piston  and  the  pressure  rise  in  the  chamber.  For 
this  shot,  these  values  were  input  to  the  model,  coming  from  the  lumped  parameter  model 
with  its  discharge  coefficient  of  0.95.  The  only  effect  of  the  low  discharge  coefficient  in  the 
model  is  that  the  cross-sectional  area  of  the  orifice  is  a  little  larger  in  the  model  than  in 
the  lumped  parameter  model.  This  is  seen  in  Figure  8,  in  that  the  injection  piston  motion 
computed  by  the  DSD/SST  model  never  catches  up  with  that  computed  by  the  lumped 
parameter  model. 

Since  the  reservoir  pressure,  the  injection  piston  motion,  and  the  discharge  coefficient 
computed  by  the  DSD/SST  model  of  the  RLPG  all  compare  reasonably  well  with  the 
expected  results,  the  rest  of  the  data  generated  by  this  model  are  considered  reasonably 
accurate.  The  goal  of  this  project  is  an  understanding  of  the  flow  phenomena  that  occur 
in  the  injector  region  of  the  RLPG.  Figures  10  and  11  show  the  Mach  number  distribution 
of  the  flow  at  6.4  ms,  8.4  ms,  10.4  ms,  and  12.4  ms.  By  6.4  ms,  the  orifice  has  not  opened 
very  much,  and  there  is  little  flow  through  the  orifice.  By  8.4  ms,  the  flow  through  the 
orifice  is  well  established,  and  a  recirculation  region  has  developed  on  the  chamber  side  of 
the  injection  piston.  The  jet  of  LP  entering  the  combustion  chamber  still  tends  to  lift  off  of 
the  nose  of  the  control  piston.  By  10.4  ms,  the  flow  through  the  orifice  has  attained  sufficient 
speed  that  the  jet  of  LP  entering  the  chamber  better  follows  the  control  piston,  although 
the  jet  still  has  detached  by  the  end  of  the  piston.  The  recirculation  region  by  the  injection 
piston  has  grown  in  strength.  At  12.4  ms,  the  orifice  is  closing,  and  the  flow  of  LP  into  the 
chamber  has  almost  stopped. 

The  jet  of  LP  lifts  away  from  the  surface  of  the  inner  piston  at  all  four  times  shown  in 
Figures  10  and  11.  The  point  at  which  the  jet  leaves  the  surface  varies  during  the  firing  cycle, 
but  much  of  the  nose  of  the  control  piston  is  in  contact  with  combustion  gases  throughout 
the  firing  cycle. 

Figure  12  show  the  Mach  number  distribution  in  the  orifice  at  the  times  previously 
mentioned.  The  highest  Mach  number  is  seen  in  the  throat  of  the  orifice,  as  expected.  The 
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Figure  10.  Shot  38:  Global  View  of  Mach  Number  Contours  at  Times  6.4  ms 
and  8.4  ms. 

recirculation  region  near  the  injection  piston  can  be  clearly  seen.  There  is  no  tendency  for 
the  flow  to  separate  from  the  pistons  or  reverse  direction  of  travel  in  the  orifice.  Either  of 
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Figure  11.  Shot  38:  Global  View  of  Mach  Number  Contours  at  Times  10.4  ms 
and  12.4  ms. 

these  could  result  in  disruption  of  LP  flow  in  the  orifice  and  possibly  cause  reversals  in 
the  gun. 
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Figure  12.  Shot  38:  Mach  Number  Contours  in  the  Orifice  at  Times  6.4  ms, 
8.4  ms,  10.4  ms,  and  12.4  ms. 

5.2  Shot  128.  Shot  128  was  fired  using  224.1  cm3  of  LP  and  an  initial  chamber  volume 
of  290  cm3  [20].  The  ignition  sequence  started  at  time  t  =  0.0.  The  simulation,  which  does 
not  include  the  ignition  sequence,  started  at  5.65  ms.  Remeshing  was  performed  every  10 


time  steps  during  this  simulation.  Figures  13  and  14  show  the  mesh  at  10.0  ms,  12.0  ms, 
14.0  ms,  and  16.0  ms.  The  meshes  have  5,602,  5,449,  8,119,  and  7,074  nodes  and  10,663, 
10,357,  15,697,  and  13,655  elements,  respectively.  As  the  firing  cycle  proceeds,  the  liquid 
reservoir  gets  smaller  in  volume.  In  the  remeshing  procedure  used  in  the  code,  the  number 
of  nodes  along  each  wall  of  the  reservoir  stays  the  same;  this  number  can  only  be  changed 
manually  when  the  code  is  restarted.  As  the  reservoir  decreases  in  volume,  the  nodes  along 
two  of  the  boundaries  get  closer  together.  Since  the  mesh  in  the  interior  of  the  reservoir  is 
based  on  the  nodes  on  the  boundary,  at  14.0  ms,  the  mesh  in  the  reservoir  is  finer  than  it 
needs  to  be.  The  remeshing  procedure  should  be  improved  to  prevent  this. 

Figure  15  shows  the  mesh  in  the  orifice  at  the  times  previously  mentioned.  The  mesh, 
especially  along  the  injection  piston,  is  sufficiently  fine  to  model  accurately  the  flow  in  the 
orifice.  Notice  that  the  orifice  opens  noticeably  more  than  in  shot  38.  This  occurs  because 
more  than  twice  the  amount  of  LP  has  to  be  injected  into  the  combustion  chamber. 

Filtered  experimental  data  was  used  as  the  input  for  the  control  piston  motion  and  the 
chamber  pressure.  The  computed  pressure  at  the  center  of  the  rear  wall  of  the  reservoir  is 
plotted  with  the  filtered  experimental  chamber  pressure  in  Figure  16.  The  reservoir  pressure 
gauge  failed  during  this  shot.  Comparing  Figure  16  to  Figure  7,  the  computed  pressure  for 
shot  128  looks  reasonable.  The  reservoir  pressure  oscillations  near  peak  pressure  are  caused 
by  a  deceleration  of  the  control  piston,  and  the  later  oscillations  are  caused  by  remeshing. 

The  ratio  between  the  reservoir  pressure  and  the  chamber  pressure  is  approximately  equal 
to  the  ratio  between  the  surface  area  on  the  chamber  side  and  reservoir  side  of  the  injection 
piston,  which  is  1.47.  Early  in  the  firing  cycle,  as  the  injection  piston  is  accelerating,  the 
pressure  ratio  should  be  a  little  lower;  later,  when  the  piston  is  decelerating,  the  ratio  should 
be  a  little  higher.  An  analysis  of  the  data  shown  in  Figure  16  verifies  that  this  is  the  case. 

The  control  piston  motion  measured  during  shot  128  is  shown  with  the  computed  injection 
piston  motion  in  Figure  17.  The  experimental  injection  piston  motion  data  from  this  shot 
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Figure  13.  Shot  128:  Global  View  of  Mesh  at  Times  10.0  ms  and  12.0  ms. 

were  not  reliable.  When  compared  with  Figure  8,  the  injection  piston  motion  from  the 
simulation  of  shot  128  is  quite  reasonable. 
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Figure  14.  Shot  128:  Global  View  of  Mesh  at  Times  14.0  ms  and  16.0  ms. 

The  discharge  coefficient  computed  during  the  simulation  of  shot  128  is  shown  in 
Figure  18.  The  flow  in  the  orifice  is  fully  developed  by  8.2  ms.  The  discharge  coefficient  drops 
below  zero  for  a  short  period  of  time  at  the  beginning  of  the  simulation.  As  the  injection  and 
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Figure  15.  Shot  128:  Mesh  in  the  Orifice  at  Times  10.0  ms,  12.0  ms,  14.0  ms, 
and  16.0  ms. 

control  pistons  separate,  there  is  a  low-pressure  region  at  the  center  of  the  orifice,  and  fluid 
from  both  the  reservoir  and  the  chamber  is  drawn  into  the  orifice.  As  the  LP  flow  develops, 
the  flow  from  the  chamber  into  the  orifice  ends  quickly.  The  discharge  coefficient  in  this  shot 
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Figure  16.  Shot  128:  DSD/SST  Prediction  of  Pressure  at  the  Back  Wall  of 
the  Reservoir  (Line).  Experimental  Pressure  in  the  Combustion 
Chamber  (Dash). 
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Figure  17.  Shot  128:  DSD/SST  Prediction  of  Motion  of  the  Injection 

Piston  (Line).  Experimental  Motion  of  the  Control  Piston  (Dash). 


is  once  again  lower  than  expected,  oscillating  between  0.8  and  0.85  vs.  the  expected  value 
of  0.95.  The  discharge  coefficient  is  increasing  toward  the  end  of  the  firing  cycle. 
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Figure  18.  Shot  128:  Discharge  Coefficient  in  the  Orifice:  DSD/SST  Prediction. 

Figures  19  and  20  show  the  Mach  number  at  10.0  ms,  12.0  ms,  14.0  ms,  and  16.0  ms. 
The  increased  opening  between  the  piston  can  be  seen.  This  leads  to  a  thicker  jet  of  LP 
entering  the  chamber.  The  jet  stays  attached  to  the  control  piston  nose  until  the  end  at 
12.0  ms  and  14.0  ms.  This  creates  a  recirculation  region  at  the  end  of  the  control  piston.  At 
the  two  other  times,  the  jet  lifts  off  of  the  piston  nose  at  different  points. 

Figure  21  shows  the  Mach  number  in  the  orifice  at  the  same  times  previously  mentioned. 
The  region  of  the  highest  Mach  number  is  in  the  throat  of  the  orifice,  as  expected.  The 
acceleration  and  deceleration  of  the  LP  entering  and  exiting  the  orifice  can  be  seen.  As  with 
shot  38,  there  is  no  tendency  of  the  flow  to  separate  from  the  piston  surfaces,  indicating  a 
well-designed  orifice. 

6.  Conclusions 

The  DSD / SST-based  finite  element  model  of  the  injection  of  the  LP  through  the  orifice 
in  the  RLPG  was  described.  The  model  was  used  to  simulate  shot  38  of  the  30-mm  RLPG. 
The  reservoir  pressure  and  injection  piston  motion  computed  by  the  model  compare  well 
with  the  values  coming  from  a  lumped  parameter  simulation  of  the  shot.  Also,  the  discharge 
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Figure  19.  Shot  128:  Global  View  of  Mach  Number  Contours  at  Times  10.0  ms 
and  12.0  ms. 

coefficient  is  reasonably  close  to  the  expected  value.  Because  of  this  agreement,  the  data 
generated  by  the  model  can  be  taken  as  at  least  qualitatively  correct.  A  total  of  three  shots 
of  the  30-mm  RLPG  were  studied. 
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Figure  20.  Shot  128:  Global  View  of  Mach  Number  Contours  at  Times  14.0  ms 
and  16.0  ms. 

Shot  38  was  a  small-charge  shot,  and  shot  128  was  a  medium-charge  shot.  The  RLPG 
performed  properly  during  both  of  these  shots.  The  flow  of  LP  through  the  orifice  tended 
to  remain  attached  to  both  pistons  through  the  orifice  and  to  the  control  piston  into  the 
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Figure  21.  Shot  128:  Mach  Number  Contours  in  the  Orifice  at  Times  10.0  ms, 
12.0  ms,  14.0  ms,  and  16.0  ms. 

combustion  chamber.  The  point  where  the  jet  lifted  off  of  the  control  piston  differed  between 
the  shots  and  varied  during  the  firing  cycle.  The  behavior  of  the  jet  in  the  chamber  is  not  to 
be  considered  very  accurate,  since  combustion  strongly  influences  the  flow  in  the  chamber. 
The  flow  in  the  orifice  behaves  well,  indicating  that  the  orifice  is  well  designed. 


This  work  provides  confidence  that  the  model  can  be  used  to  study  the  flow  of  LP  from  the 
reservoir  and  through  the  orifice.  However,  since  it  neglects  combustion,  the  model  cannot 
predict  phenomena  that  occur  in  the  chamber.  A  current  project  involves  the  coupling  of 
the  reservoir  model  with  a  model  of  the  chamber,  which  includes  combustion.  The  resulting 
model  will  be  used  to  study  the  entire  firing  cycle  of  the  RLPG. 
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